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We study the normal and the superconducting properties in noncentrosymmetric heavy fermion 
superconductors CeRhSi3 and CeIrSi3. For the normal state, we show that experimentally observed 
linear temperature dependence of the resistivity is understood through the antiferromagnetic spin 
fluctuations near the quantum critical point (QCP) in three dimensions. For the superconducting 
state, we derive a general formula to calculate the upper critical field H C 2, with which we can 
treat the Pauli and the orbital depairing effect on an equal footing. The strong coupling effect for 
general electronic structures is also taken into account. We show that the experimentally observed 
features in H C 2 \\ z, the huge value up to 30(T), the downward curvatures, and the strong pressure 
dependence, are naturally understood as an interplay of the Rashba spin-orbit interaction due to the 
lack of inversion symmetry and the spin fluctuations near the QCP. The large anisotropy between 
H C 2 || z and H C 2 _L z is explained in terms of the spin-orbit interaction. Furthermore, a possible 
realization of the Fulde-Ferrell- Larkin-Ovchinnikov state for H _L z is studied. We also examine 
effects of spin-flip scattering processes in the pairing interaction and those of the applied magnetic 
field on the spin fluctuations. We find that the above mentioned results are robust against these 
effects. The consistency of our results strongly supports the scenario that the superconductivity in 
CeRhSi3 and CelrSi3 is mediated by the spin fluctuations near the QCP. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

In noncentrosymmetric heavy fermion superconductors, in addition to strong electron correlation, there exists an- 
other key property, the anisotropic spin-orbit (SO) interaction due to the lack of inversion symmetry. The anisotropic 
SO interaction plays important roles both in the normal and the superconducting state, and is expected to lead to 
many interesting phenomenal— For such phenomena, electron correlation is quite important, because it can largely 
enhance the effect of the SO interaction. The interplay of the anisotropic SO interaction and electron correlation is 
truly an unique nature in noncentrosymmetric heavy fermion compounds.— ~— In particular, such an interplay in the 
superconducting state has been attracting particular interest. In this context, especially, CeRhSia^r— and CeIrSi j 16 i 17 
are promising candidates for the interplay, because they are considered to be located near the antiferromagnetic (AF) 
QCPs around which strong correlations through the spin fluctuations are essential. 

CcRhSi3 and CeIrSi3 are AF metals at ambient pressure, and begin to exhibit superconductivity at some critical 
pressures P c where the Neel temperatures seem to be suppressed to absolute zero. According to the neutron experi- 
ments for CeRhSi3, the AF ordering vector is Q = (±0.437r, 0, 0.57r), (0, ±0.437r, 0.57t) and the nature of the AF order 
is spin density wave-like,^- This is different from CePtaSi in which the AF seems to have localized nature and the su- 
perconductivity coexists with it even at zero applied pressured In NMR experiments in CeIrSi3, 1/Ti oc Tj (T + 9) 1 / 2 
is observed near the critical pressure, which is a characteristic behavior of the systems with 3-dimensional (3D) AF 
spin fluctuations <22r— In addition, the resistivity p in both CeRhSi3 and CeIrSi3 above the superconducting transition 
temperatures T c in some pressure regions near the QCP shows the anomalous T-linear dependence which is different 
from p ~ T 2 in canonical Fermi liquidsi 15 ' 16 

The QCP related phenomena are observed also in the superconducting state. The large jump in the heat capacity 
at T c in CeIrSi3 can be attributed to the strong coupling effect due to the spin fluctuations^ It has strong pressure 
dependence and is largely enhanced near P c . The most striking phenomena which would be related to the quantum 
criticality in CeRhSi3 and CeIrSi3 appear in the behaviors of the upper critical fields H C 2 when the applied magnetic 
field is parallel to z-axis i 24 ' 25 The remarkable features of the experimental results are as follows, (i) As the pressure 
approaches a critical value, H C 2 exhibits extremely high value which exceeds the orbital limit as well as the Pauli limit 
estimated by the conventional BCS theory. The observed ~ 30(T) is the highest value among the heavy fermion 
superconductors ever discovered, although T c is merely T c ~ 1(K). (ii) Hc2 curves have downward curvatures and 
the increase is accelerated as the temperature is decreased, making a sharp contrast to any other superconductors in 
which the increase of H C 2 becomes slower as T is decreased, (iii) H C 2 increases very rapidly as the pressure approaches 
the critical value, while the pressure dependence of T c is moderate. These characteristic features strongly suggest that 
there exists a deep connection between the superconductivity and the magnetic quantum criticality. In the previous 
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study, the present authors have shown that these experimental results are well explained as an interplay of the Rashba 
SO interaction due to the lack of inversion symmetry and the spin fluctuations near the QCP^ 

On the other hand, H C 2 for in-plane fields differs from H C 2 || z in some important features. H C 2 _L z is merely 
less than 10(T) and its pressure dependence is moderate, and the H C 2 curves exhibit usual upward curvatures j 24 ' 25 
This anisotropy in H C 2 would be related to the Rashba SO interaction, since the Fermi surface is asymmetrically 
distorted by the in-plane field and the Pauli depairing effect plays essential roles. By contrast, the renormalization of 
the quasiparticle velocity by the spin fluctuations is almost isotropic, resulting in the enhanced orbital limiting field 
in all directions of the applied field. Another interesting phenomenon in the noncentrosymmetric superconductors in 
applied magnetic fields is the helical vortex phase which has been discussed theoretically. —— In a helical vortex 
phase, the superconducting gap function is modulated in real space, A(J?) ~ cxp(zQ • R)A with the modulation 
vector Q ~ {a/ 'ef)pbH /vp, where a,vp,Sp and p,s are the strength of the SO interaction, the Fermi velocity, the 
Fermi energy, and the Bohr magneton, respectively. For 3D Rashba superconductors for H _L z, however, it is pointed 
out that this phase modulation is just a translational shift of the vortex lattice and has no physical importance^ - — 
Several authors also have discussed a spatially modulated superconducting state under magnetic fields with a large 
Q ~ hbH/vf which is continuously connected from Q ~ (a/eF )/j,bH/vf - 2S]29]31]33 This large Q state corresponds 
to the Fulde-Ferrell-Ovchinnikov-Larkin (FFLO) state i 35 ' 36 The stability of the modulating superconducting state 
with the large Q depends on the relative strength of the orbital depairing effect to the Pauli depairing effect in the 
compounds. 

In this paper, we study the normal and the superconducting properties in noncentrosymmetric superconductors 
CeRhSi3 and CeIrSi3. We examine the anomalous T-linear dependence of the resistivity in the normal state. In the 
previous studies ; 21 ' 22 at very low temperatures, p ~ T 3 / 2 is predicted for 3D AF spin fluctuations. The temperature 
dependence of the resistivity for 3D AF spin fluctuations have been studied in detail by several authors 3 ^- 41 - We, here, 
show that p ~ T in CeRhSi3 and CeIrSi3 is actually due to the 3D AF spin fluctuations. For the superconducting state, 
the upper critical fields both for H \\ z and H _L z are investigated. For the calculation of H C 2, we derive a general 
formula which enables us to treat the Pauli and the orbital depairing effects on an equal footing. We can also take into 
account the strong coupling effect for a given electronic structure. We calculate H C 2 on the basis of the scenario that 
the superconductivity in CeRhSi3 and CeIrSi3 is mediated by the spin fluctuations, and show that the experimental 
features are well explained as an interplay of the spin fluctuations and the Rashba SO interaction. Although the 
formula is applicable for general models, we use a phenomenological model to calculate H C 2 and neglect the following 
two points in the model. One is the scattering processes in the pairing interaction in which spins of quasiparticles 
are flipped by the Rashba SO interaction. It is pointed out that such processes can enhance the admixture of the 
singlet and the triplet superconductivity ) 37 ' 38 and the strength of the admixture affects H C 2 _L z We show that the 
admixture is still small in CcRhSi3 and CeIrSi3 even if we include the spin-flip scattering processes. The other point 
is the applied field dependence of the spin fluctuations. Because the applied field is so large in CeRhSi3 and CeIrSi3 
especially for H \\ i that one might think that the spin fluctuations are suppressed and they cannot contribute to the 
enhancement of H c2 . We show that the spin fluctuations are robust against the applied field H up to the strength of 
the Rashba SO interaction, (IbH < a, because the Rashba SO interaction tends to fix the directions of spins on the 
Fermi surface and it competes with the Zeeman effect. The consistency of our results with the experiments strongly 
supports the scenario that the superconductivity in CeRhSi3 and CeIrSi3 is mediated by the spin fluctuations near 
the AF QCP. 

This paper is organized as follows. In Sec|TTl we study the experimentally observed T-linear dependence of the 
resistivity. In Sec lIIIi a general formula for the calculation of H C 2 is derived from the Eliashberg equation. The 
characteristic features of H C 2 in CeRhSi3 and CeIrSi3 are well explained with the use of the formula in Sec lIVl We 
discuss, in Sec|Vl the spin-flip scattering processes in the pairing interaction and the magnetic field dependence of 
the spin fluctuations which are not included in the approximation used for the computation of H C 2- The summary is 
given in Sec lVIl 

II. RESISTIVITY IN NORMAL STATE 

In this section, we discuss the temperature dependence of the resistivity near the AF QCP in CeRhSi3 and CeIrSi3. 
In CeIrSi3, NMR 1/Ti behaves as 1/Ti cx X ' j\jT + 6 in some pressure regions, which means that the character of 
the spin fluctuations is 3D antiferromagnetic^r— In noncentrosymmetric systems, however, spin fluctuations are not 
isotropic due to the anisotropic spin-orbit interaction. The anisotropy in the noninteracting susceptibility \° is of 
the order of a/ Sf <C 1, where a is the strength of the spin-orbit interaction and Ef is the Fermi energy. Actually, in 
Sec lVBl we show that the anisotropy among Xxx,Xyy an d Xzz is very small within the random phase approximation. 
Therefore, we can neglect the anisotropy of the spin fluctuations for the discussion of the resistivity and H C 2 in 
CeRhSi 3 and CeIrSi 3 . 
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For the systems with 3D AF spin fluctuations, the resistivities are expected to be p ~ T 3 / 2 according to the 
previous studiesi 21 ' 22 In CeRhSi 3 and CeIrSi 3 , however, the temperature dependence is p ~ T near the AF QCPi 15 i 16 
The resistivity due to the 3D AF spin fluctuations was discussed by several authors, and p ~ T behavior was found 
for the weakly disordered system s 39 ' 40 and the clean systems4i Here, we show that the T-linear resistivity in CeRhSi3 
and CeIrSi3 is naturally understood in terms of the 3D AF spin fluctuations, and this behavior has basically nothing 
to do with the lack of inversion symmetry. 

CeRhSi3 and CeIrSi3 are heavy fermion systems with Kondo temperature Xk ^50-100(K) which is much higher than 
the superconducting transition temperature T c ~ l(K)i 13 i 16 i 42 Therefore, to study the properties at T =1-10(K), we 
consider the low energy quasiparticles mainly formed by /-electrons through the hybridizations with the conduction 
electrons. We use the following single band model for the low energy quasiparticles with the asymmetric spin-orbit 
interaction 

S = So + Ssf, (1) 
S = c\\-iu n + £o(fc)]c fe + 4a£o(fc, -H") ■ crc k , (2) 

k k 

Ssf = ~^^x{q)° aa < ■<?/}[)' cl +qa c ka .cl,^ q pc k .fi,i (3) 

kk' q 

where cfi is the annihilation (creation) operator of the Kramers doublet of the Tj state. Ssf is introduced phe- 
nomenologically and represents the interaction between the quasiparticles by the strong spin fluctuations near the AF 
QCP. Since CeRhSi3 and CeIrSi3 have body-centered tetragonal lattice structures with lattice spacing 1:1: 2 ) 14 i 17 
the dispersion relation £o(fc) and the Rashba type SO interaction are approximated by 

£o(k) = — 2ti (cos k x a + cos k y a) + 4*2 cos k x a cos k y a 

—8*3 cos(k x a/2) cos(k y a/2) cos k z a — 2*4 cos 2k z a — p, (4) 
C (k,H) = (sinkyci, - s'mk x a - fi B H v /a, -p B H z /a), (5) 

where a is the lattice constant and p is the chemical potential. Although H — in this section, we include the Zeeman 
effect in the action for the later discussion. We fix the parameters as (*i, *2, *3, *4, n, a) = (1.0, 0.475, 0.3, 0.0, 1.0, 0.5) 
by taking t\ as the energy unit. The Fermi surface determined by these parameters is in qualitative agreement with the 
band calculation and can reproduce the peak structures of the momentum-dependent susceptibility observed by the 
neutron scattering experiment o 18 ' 43 ' 44 . Since we consider /-electron systems, we assume that the above parameters 
include effects of the mass renormalization due to local spin correlations with typical energy scale Tk ~ 50-100 (K), 
the Kondo temperature; i.e. t\ ~ 50-100 (K). 

The interactions are phenomenologically introduced through the renormalized susceptibility X^) ) 21 ' 22 ' 26 ' 45 ' 46 

X(iv n ,q) = T ^2 ( 6 ) 



where %oj <Zo an< i T are respectively the susceptibility, the length scale and the energy scale of spin fluctuations without 
strong correlations. These quantities are renormalized through the coherence length £(T) as the system approaches the 
QCP. The critical exponent of £ is the mean field value 1/2 and 9 is considered to decrease monotonically as the applied 
pressure approaches the critical value for the AF orderi 21 i 22 The temperature dependence of £ is also consistent with the 
recent NMR experiment for CeIrSi 3 .^°- The ordering vectors are Q l = (±0.43?r, 0, 0.5?r)/a, Q 2 = (0, ±0.43?r, 0.5?r)/a 
according to the neutron scattering experiments for CeRhSi3.— In this study, we fix the parameters in xfan, <z) as 
go = a~\ r = 3.6*i and £ = 0.45a. The value T = 3.6*i is of the same order as the Fermi energy without the effect 
of the spin fluctuations. £ is determined so that the maximum of £ would be £ max < 10a at the lowest temperature 
in this study, which is a reasonable value for the AF spin fluctuations. We note that the coupling constant g should 
be regarded as an effective one renormalized by the vertex corrections . 47 ' 48 
The Green's function isJ^r— 



G a p(k) = lr a0 (k)G T (k), (8) 

r=±l 

W(*0 = I ( 1 + T ^( fc ) • °) . ( 9 ) 



iu n - s T [k) - S (k) 
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where e T (k) = e (k) + Ta\\C(k)\\, C{k) = C (k) + H(k)/a, C(k) = C(k)/\\C(k)\\ and ||£(fc)|| = yEti^Wl 2 - 
Selfenergy is introduced as E = (E n + E u )/2, S x = (E it + E n )/2, E y = (£ it - S n )/2i, and E z = (S tt - S u )/2. 
Up to the first order in g 2 Xo, an d S are expressed as 

^o(k) = ^Y,9 2 x(k~k')[G° n (k') + G a u (k% (11) 

k' 

= ^Ey^ fc - fc ')[-G; t (fc')-G?,(fc')], (12) 

k' 

^y(k) = ^El^-^K^ + ^l (13) 

fc' 

= ^E^-^K^ + GJ;^)] (14) 

k' 

where G°^ is the non- interacting Green's function. We have neglected the constant terms in E^. In the above 
expression of E M , the most dominant term is Eo, and ^x,y is smaller than Eo by the factor cx/sf <C 1, where Ef is 
the Fermi energy. For Rashba superconductors, E z = without magnetic field. 
The conductivity is calculated from the Kubo formula 

o> = Iim -Imi^o;), (15) 
K^(iuj n ) = I dTe tUnT (TJ^T)J v (0)}, (16) 



o 



where the current operator J„ is defined as 



J» = ey^c^fc^Cfc, (17) 

k 

Vk„ = V^[e (k)+aC Q (k)-cr}. (18) 

After the analytic continuation, G A G R , has the dominant contribution to the conductivity. Among the four terms 
{G A G R ,} T . T ' = ±i, for sufficiently large a, the terms G A G R and G^G+ have no singularity with respect to the quasi- 
particle damping rate. Therefore, we can neglect them, and the resulting expression for a xx is 

£tr[« te &«jj&] (-%(e k r)j (19) 



e 2 



k. 



7 fer = -[ImS^(0, fe) + rRe£ fl (0, fe) • lm£ fl (0, k)} (20) 

where Re£ R = ReC R /\\ReC R \\, l [ k R ' A) = li R ' A) (s kr ,k) and e kT = e T {k) + ReE (0,fc). Here we have neglected the 
vertex corrections which are necessary for the current conservation law.— ~— This is because, for the resistivity, the 
scatterings by the AF spin fluctuations have large momentum transfers, and therefore, the back-flow included in the 
vertex corrections does not affect the temperature dependence of the resistivity^ As mentioned above, X is much 
smaller than Eo in amplitude, and from Eqs. (ITT|) ^ dT4l . the temperature dependence of Eo and that of X are the same. 
Therefore, hereafter, we neglect S and take into account only Eo in this study. The resistivity in noncentrosymmetric 
systems is almost the same as that in usual centrosymmetric systems. This is different from the situations for the 
anomalous Hall effect, the magnetoelectric effect and so on for which the Rashba SO interaction plays important 
roles^iSria 

Before moving to the numerical calculation, we show a simple analytical result for ImEo which determines the 
qualitative behavior of a^ u . Since Eq. (|ll[) is basically the same as the selfenergy E cen in the usual centrosymmetric 
systems, we consider E cen for brevity. The selfenergy at the hot spots for a sufficiently clean system with the impurity 
damping "fi mp -C vf/S, is calculated as 

27 r° th 2T " tanh 2tJ Imxi?(e ' 9) ImG -n(e> *f - <z) 

^ XqT s (itT) 2 oj? 

1^ 7777, _j_ Tr r j~ ] i ecu ( ' k F - q) 

~ e°rin[l + 7rT/2r s ], (21) 
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where G 0R is the retarded Green's function including the impurity damping 7i mp and \ R is the retarded susceptibility 
obtained from the analytic continuation of Eq. © . In the above calculation, the dispersion at the hot spots has been 
expanded as eo(kp — q) = \7eo(kp — Q) ■ (Q — q) for q~ Q, because £o(kp — Q) = —eo(kp) = is satisfied at the 
hot spots. Here, we have used the approximation H(z) = J c?e[coth(e/2T) — tanh(e/2T)]e/ {uj 2 + e 2 ) = l/z — 2[tp(z + 
1) — tp(z + 1/2)] ~ (ttT) 2 /[uj q (uj q + ttT/2)], where tp is digamma function and z = uj q /(2TrT)i 51 ^ 3 This approximate 
form becomes exact both for z — > and z — > oo. %Q,r s and uj q are defined as \Q = Xq(£,<1o) 2 ,F s = To(t;qo)~ 2 and 
— r s [l + £, 2 (q — Q) 2 ], respectively. Therefore, we have p ~ — ImS^ ~ T when the hot spots are dominant for the 
conductivity. This is a general behavior for the clean 3D systems with the 3D AF spin fluctuations. 

In the numerical calculation, we neglect the real part of the selfenergy which changes the shape of the Fermi 
surface, because such an effect is non-perturbative. We regard e T (fe) as the dispersion that includes ReEo- We show 
the numerical results for p xx — l/a xx by using Eq.fTTj) and Eq. (fT9|) for clean limit. As shown in FigJU for sufficiently 




FIG. 1: Resistivity v.s. temperature for several 6. From the top to the bottom, 9/ti = 0.002, 0.05, 0.2. 



small 0, the resistivity is proportional to T in a wide range of temperature where the hot spots are thermally blurred 
and dominant for the conductivity. In a very low temperature region where such blurring is suppressed, p is dominated 
by the electrons in the cold spots. For large 9, the canonical Fermi liquid behavior p ~ T 2 can be seen. The calculated 
p well explains the experimentally observed features of the resistivity in CeRhSi3 and CeIrSi3. Therefore, we conclude 
that the observed p ~ T above T c is due to the AF spin fluctuations. 

We put a remark on the impurity effect i 39 ' 40 If the impurity scattering is sufficiently strong, the anisotropic scatter- 
ings by the spin fluctuations are smeared, which weakens the singularity. We, here, simply estimate the selfenergy by 
the spin fluctuations in the presence of the impurities for centrosymmetric systems analytically. For the system 
with the strong impurity effect which smears the anisotropy by the AF spin fluctuations, we evaluate the selfenergy 
E^ n averaged on the Fermi surface 

(ImI] ccn (0,fe)) FS = EfcImG o« (0jfc) 

v x Q r,(7rr) 2 aim X °£(o,g) 

^ ^ La q (w q + -KT/2) doj 

~ r^Vl+wT/Zr.-l), (22) 

where we have defined Xconl?) = ~ (T/N) J2k ^ccn(k)Gc Cn (fc + q). Here, we have assumed that its T, g-dependence 
is moderate and it does not contribute to the selfenergy. We obtain p ~ T 3 / 2 for £~ 2 ~ T in dirty systems. Note 
that, in the case of - T 3 / 2 , we again have £- x T ^1 + nT/2T s - lj ~ T i / i T l T~ 1 / 4 ' = T 3 / 2 for sufficiently low 

temperatures. Thus, the resistivity in the dirty systems with 3D AF spin fluctuations is p ~ T 3 / 2 both for £~ 2 ~ T 
and £~ 2 ~ T 3 / 2 in agreement with the previous studies ! 21 ' 22 



III. ELIASHBERG EQUATION IN MAGNETIC FIELD 



A. exact formula within semiclassical approximation 



In this section, we derive a formula for the calculation of H C 2 from the linearized Eliashberg equation in real 
space. The derivation is based on the semiclassical approximation which is legitimate for the systems with kpla 3> 1, 
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where kp is the Fermi wave number and Ih — l/yjefff is the magnetic length. This condition is satisfied for many 
superconductors including heavy fermion superconductors, and therefore, the resulting equation for H C 2 is applicable 
for a number of compounds. Our formula is a generalization of the previous studies,— ~— and can be extended easily 
to more complicated models although we use a single band model in this section. 

To derive the formula for the calculation of H C 2, we use the linearized Eliashberg equation in real space with the 
vector potential A which gives a uniform magnetic field, 

A aa ,(iuj n7 x,x';A) = -— ^ V aa >^^ (w ro , x, x'; iu m , y, y'; A) 

iu m yy 

x Gf3~,(iu m ,y, z; A)A lY (iu m , z, z'; A)Gp,^,(-iu m , y' , z'\ A), (23) 

ZZ' 

where represents the summation over all lattice sites, and the spin indices are summed over. G, A and V are, 
respectively, the normal Green's function, the gap function and the pairing interaction. Note that, if A is fully 
taken into account in the above equation, the resulting equation is gauge invariant under the gauge transformation 
A(x) —> A(x) + Wf(x) and i/j(t,x) — > exp[ief(x)]if>(r, x) where ip ls the field operator of the electrons. By this 
transformation, each factor in the equation acquires the additional phases as, 

G(ioj n ,x,x') -> exp[ie(f(x) - f(x'))]G(iu n ,x,x'), (24) 
A(iu n ,x,x') -> cxp[ie(f(x) + f(x'))]A(iuJn,x,x'), (25) 
V(iuj n ,x,x';iuj m7 y,y') -> exp[ie (f(x) + f{x') - f(y) - f(y'))]V(iui n , x, x'; iuj m , y, y'). (26) 

In this study, however, we use the semiclassical approximation in which we do not explicitly include the effect 
of the vector potential on the pairing interaction V, because the vector potential in V is not responsible for the 
Landau quantization of the gap function which is the most important phenomenon of the orbital effect in type-II 
superconductors. By contrast, the lack of translational invariance in G and A in the presence of the applied vector 
potential A is related to the Landau quantization. Within the semiclassical approximation, the normal Green's 
function is 

G(iuj n ,x,y;A) = e^ x '^G{iu n , x - y; A = 0), (27) 
<p(x,y) = e / A{s)ds. (28) 



We can easily perform the integral along the straight line s(t) = y + t(x — y),0 < t < 1, using the relation 
A(ax + by) = aA(x) + bA(x) which holds for any A giving a uniform magnetic field H, and obtain 

<p(x,y) = eA(^-^-).(x-y). (29) 



Although the linearized Eliashberg equation is no longer invariant under the gauge transformation defined above 
within this approximation, it is still gauge invariant under a gauge transformation which involves only the center of 
mass coordinate of the Cooper pairs. 

Next, we proceed to rewrite Eq. (|23[) in fc-space. The pairing interaction V should be decomposed into two parts as, 

V(%un,x, x'; iu m , y, y') = V re \ioj n ,x - x'; iu, m , y - y') x V ccn (^-^; V±jC\ (30) 



where V rcl and V ccn are the interactions in the relative coordinate and the center of mass coordinate. Here, we take 
V cen to be dimensionless. It is convenient to introduce the following variables 

_ x + x' 
It - 



2 ' 

Y = y-z, Y' = y'-z', (31) 

. z + z' , , 
R' = — — , r' = z-z'. 
2 

In this coordinate, the phase factor exp(icp(y, z) + iip(y', z')) which arises from G(iu) m , y, z)G(—iu) m , y' , z') in Eq. ([2^|) 
becomes 

<p(y, z) + tp(y', z') = eA(R')(Y + Y') + eA(r')(Y-Y') + eA(Y)Y + eA(Y')Y' 
~ 2eA(R') Y+ n Y + eA(Y)Y + eA(Y')Y' . 
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In the second equality, the neglected term eA(r')(Y — Y') is much smaller than the first term, since for the dominant 
scattering processes, \\Y — Y'\\, \\r'\\ <§C \\Y + Y' || are satisfied in the systems with short range pairing interaction. 
The first term represents the phase which the Cooper pair with center of mass R' acquires. We perform the Fourier 

transformation of Eq. (|23f and assume V^™^, (^R; R' + ^ = <5r„r/+(y+y')/27 then we obtain 



kk' PP' YY 



Y + Y' 



x Gp 7 (iuj m ,p)Gp>y (-iu m , p' ) A 77 / (iuj m ,k',R 
x exp i (kr + (-k' + p)Y + (k' + p')Y') 

x expi (eA(Y)Y + eA(Y')Y' + ca(r - Y + Y \ (y + Y'\ 



The phase factor including A is rewritten as 

/ f Y + Y'\ Y + Y'\ ( Y + Y' 
exp % ( eA(Y)Y + eA(Y')Y' + 2eA ( R - J - J A 77 > ( iu m , k',R - 



Y + Y' 

expi ( n(H) ) A 7 y (iuj rn ,k\R 



where Tl(R) = -iV -2eA(R), 1 = -e[(Y + Y')/2} 2 V R A(R) and 9 2 = (e/2)A(Y - Y')(Y - Y'). 6 2 is proportional 
to (y — Y'), and therefore, negligible. 9\ is also small compared with U(Y + Y')/2, because 9\ ~ lKkpln) 2 while 
U(Y + Y')/2 ~ Neglecting 6>i and 02, we end up with the linearized Eliashberg equation in fc-space in the 

presence of the vector potential A, 

A aa ,(k,R) = -^^V r aa /, W /(fc,fc / )G fl7 (fc' + n/2)G fl v(-A;' + n/2)A 7 y(A;',il), (32) 

where k = (ioJ n , k) and II = (0, II), and we have written V rcl as V for simplicity. This is a well-known form of the 
Eliahsberg equation and similar expressions are often used for the discussion of H c2 in superconductors. As mentioned 
before, if we define a semiclassical gauge transformation which involves only R as 

A(/c, R) -> exp[z2e/(H)] A(fc, R), (33) 

this equation is gauge invariant, because, for 0[A{R)] = G(k + U/2)G(-k + U/2), e~ i2e fO[A + V/]e l2e/ = 0[A] is 
satisfied. The relative coordinate is not involved in the gauge transformation in the semiclassical approximation, and 
V rel (k, k') does not change under the transformation. 

We, next, proceed to rewrite the above Eliashberg equation to perform numerical calculations. In the present study, 
we denote the coordinate as (Ri, R 2 , R 3 ) = (R x ,R y ,R z ) for the perpendicular field and (R\, R 2 , R3) = (R x , R z , R y ) 
for the in-plane field. With this notation, the gap function for H = He 3 is expanded by the Landau functions, 

00 

A aa ,(k,R) = A QQ 'n(fc)0AQn(fi), (34) 

<t>KQn{R) = e l Q ARA MRi,R2), (35) 

where {</>„} are the usual Landau functions, R A = {A 1 ^ 2 R 1 ,A' 1 / 2 R 2 ,R 3 ), and Q A = (A" 1 / 2 ^, K 1/2 Q 2 , Q 3 ). The 
parameters Q and A represent, respectively, the modulation of the gap function and the anisotropy of the vortex 
lattice in the R\R 2 plane, and both of them are optimized to give the largest H c2 . We introduce the operator 

n AQ = (-iV A - 2eA(R A )j - Q A with V A = d/dR A . The Landau functions {</>aq„} satisfy the following relations, 

U AQ (t) A Q n = Vfl + 10AQn+l, (36) 
H- Q <f>AQn = Vn<f>AQn-l (37) 



where II± Q = l -f (ll AQ qp ill 



AQ 
2 
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By taking an inner product of Eq. ([32|) . we obtain 

A aa / n (fc) = -^^^^(fc^O^/r/ST^O^^vC-^O^rrwCfcOATY^'.^)' ( 38 ) 
k' tt' 

oo 

Q Tl T 2 n in2 {k,H) = ^ (4>AQ ni I G Tl (fc + II/2) 1 0AQm) {4>KQm \ G T2 (- fc + II/2) ] (^AQ W2 ) , (39) 
m=0 

where the completeness relation ^ m |0AQm)(</>AQm| = 1 is used. Here, we have neglected the II operators in l T because 
they only lead to the terms with positive powers of eH, i.e., l T (k + II) ~ l T (k) + VZ T (fc)II = 0(1) + Ofi/kpln), while 
Q is proportional to (|e|H) -1 / 2 describing the non-perturbative effect of the formation of the vortex lattice. 

Within the semiclassical approximation, Eqs. (|3"5)l and ([3TR) are exact. For numerical calculations, however, we 
need a cut off in the summation J2m=o which should be large enough for the calculated results to be reliable. It is 
hard to solve Eq.([38f and (|39|) with such a large cut off. So, we introduce an alternative formula for the numerical 
calculation of H C 2 in the next section. 



B. alternative formula for numerical calculation 



As mentioned at the end of the previous section, it is difficult to solve the exact formula Eqs. (|3"5f and (|39l) numeri- 
cally. Then, we approximate them by an alternative equation. Instead of Eq. (|32[) . we introduce a modified Eliashberg 
equation, 

A aa ,(k, R) = -^Lj2 Vaa-Wik, k')[G^{k' + n)Gp,y(-k r ) + G M {k')G ^>{-k' + II)]A 7 y (*', R). (40) 
k' 

This equation is rewritten as 

A aa , n (k,R) = -^^y QQ ^^(fc,fcO^/r/3 7 (fc^^)/r'/Jv(- fc ^^)^rw(fc',ff)A 7y (fc',i^), (41) 

k' tt' 

g Tl T 2ni n 2 (k,H) = i[(0 AQ „ 1 |G ri (A: + n)|0AQn 2 )G r2 (-fc)+G Tl (fc)(0AQn 1 |G T2 (-fc + n)| ( /.AQ„ 2 )]. (42) 

This is the alternative formula for numerical calculations, and does not need an infinite summation like in 
Eq. (|3"9")l . Equation (|3"2"|) is not exactly equivalent to Eq. (j4"0|) when V and A are fc-dependent as in unconventional 
superconductors. However, we have confirmed that the two different formulae give the qualitatively same results for 
H C 2, and the quantitative difference is small. Therefore, hereafter, we use Eqs. (j4Tj) and (142]) . 

With the use of the relation A = J 00 dte~ at for Re(a) > 0, the matrix elements are calculated as 

min{ni,n2} 



AQni \ G Tl (k + n)|0AQn 2 ) £ 7 ^ Wyj f tf-T^-"' / ^-^-1,^+^-21 



mm{ni,n 2 } , ; - y , . "i+^-at+i 



mirH 



,_ n (nx - 0!(n 2 - 1+ W 



^AQnJG^^fc + n)!^^) = is T2 ( ni - l)\(n 2 ~ Z)!?i! C ^+ Vaz") ^ni+n 2 -2i(22) (44) 



where 



2V / \ 

f n(z) = El h^w. (45) 



n=0 x 

/• oo 



f n (z) = e z / dte-' t n . (46) 
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The variables a, 6, c and z are given by, 

ai =A Tl (k;H) = l H 2 [v^ i± (k;H)} 2 

h = s Tl fw n (fc; IT) + i[s Tl (k; H) + (k; H) ■ Q A 



zi = bx/y/2al, 

ci± = C Tl± (k;H) = -is^v^^H), (47) 
s ri = sgn(u; Tl (fc;if)), 
v A ± {k;H)=v A l {k;H)±iv A 2 {k;H), 



and 

a 2 = A T2 (k;~H), 

b 2 = s T2 (oj T2 (k; -iT) + i[-£ T2 (fc; -JT) + <(fc; -if) • Q A ]) , 

z 2 = fc/v^, (48) 
C2± = C T2 ±(k; -H), 
Sr 2 = sgn(cj T2 (fc;-F)), 

where w T (fc; H) = u n - ImS (fc, ff), e T (fc, #) = £ T (fc, H) + ReS (fc, #), and v A (k, H) = V A e r (fc, if). 

A convenient expression of f n is obtained through the recurrence formula which is directly derived from Eq. (|46[) , 



fo(z) = ^e^erfc( 2 ), (49) 

h(z) = \, (50) 

n — 1 1 

/«(*) - ^-/n-aW - 2* = ("> 2 )- ( 51 ) 



The solution is 



fc=i v ; 

where g„ = §(n : even), ^^(ra : odd) and r n — 0(n : even), l(n : odd). 

With the expression Eqs. (|43l) and (|44[) . the numerical calculation of Eq. fiTj) is straightforward. Similar expression 
for Q can be obtained in the same way and we can also solve Eq. (|38p numerically. As mentioned above, Eq. (|38[) and 
Eq. ([4"Tj) give the qualitatively same results and the quantitative difference is small. 

The important point is that Eq. (|4"Tj) allows us to calculate H C 2 for general lattice models with arbitrary Fermi 
surfaces, taking into account both the orbital and the Pauli depairing effect on an equal footing. In the Gintzburg- 
Landau approach, the relative strength of the orbital and the Pauli depairing effect is characterized by the Maki 
parameter a m = \^2H or b/Hp, where -f/orb and Hp are the orbital and the Pauli limiting field, respectively. In our 
formulation, however, we do not need such a parameter which is difficult to be determined experimentally. The 
parameter corresponding to a m in this study is an effective mass of the quasiparticle for the cyclotron motion 

Welt = " 9, (53) 

where t\ is the energy unit of the lattice model and a is the length unit, i.e., the lattice constant. Writing Ijj = Ihcl 
and (J-bH = ht\ with dimensionless variables Ir and h, we have a simple identity, 

= M$JL. (54) 
fi B m ef [ 

A large effective mass corresponds to a slow velocity of the quasiparticles for the cyclotron motion leading to a 
suppression of the orbital depairing effect. m e g can be determined reasonably, while evaluating ajj from experiments 
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is rather difficult because i/ rb and Hp are not directly observed, especially for the strong coupling superconductors. 
The lattice constant is determined experimentally, and we fix a = 4.0(A) in this study, which is consistent with the 
experiments . 14 ' 17 We can also determine the value of t\ in a reasonable way. By solving the Eliashberg equation (j40|) 
at H = 0, we obtain T c = T c t\ in the unit of t\. Then, comparing it with the experimentally observed transition 
temperature T£ xp , we have 

T^exp 

h = -jf- (K). (55) 

In this way, the parameters of the model are evaluated. However, the choice of all the parameters is not unique 
and there remains some ambiguity especially for the strength of the interaction on which T c — T c t\ largely depends. 
Therefore, we change the value of the strength of the interaction depending on the choice of the magnitude of t\ to 
make T c consistent with the observed values. For the calculation of H c2 in CeRhSi3 and CeIrSi3, we use two values 
of ti and compare the results. 

In addition to the treatment of the Pauli and the orbital depairing effect, the strong coupling effect can be included 
naturally in Eq. (|4ip . Once we calculate the pairing interaction V and the selfenergy E for a given Hamiltonian, they 
are directly incorporated into the Eliashberg equation (l4"Tj) . This feature is essentially important for the study of H c2 
in CeRhSi3 and CeIrSi3, because it is considered that they are located near the AF QCPs and the quasiparticles 
interact with each other through the strong spin fluctuations. 



IV. CALCULATION OF UPPER CRITICAL FIELD 



In this section, we show the numerical results calculated from Eq. (|4"Tj) with Eg. (14"21) . We solve the Eliashberg 
equation both for H = and H ^ 0. For the latter, we study the two cases: H = (0,0, H) and H = (0,H,0). In 
the case of H \\ z, the Pauli depairing effect is strongly suppressed by the anisotropic spin-orbit interaction, and H c2 
is determined by the orbital limiting field i/orb^ 1 11,12 On the other hand, for H _L z, the Pauli depairing effect is 
significant because of the anisotropic distortion of the Fermi surface due to the Rashba SO interaction, and H c2 is 
mainly determined by the Pauli limiting field Hp. 

In this section, we use the action Eq.(JT]) to calculate H c2 . The selfenergy Eo has the real and the imaginary part 
which have different effects, respectively. ReEo only gives the deformation of the Fermi surface, and, as in Sec|TTl it 
is reasonable to consider that e T (k) already includes the shift due to ReEo and to replace e T (k) + ReEo — > £ T (k)- On 
the other hand, ImEo(iw n , fe) gives two important effects for the quasiparticles around the Fermi level. One is the 
damping factor 7 = — ImE^ and the other is the mass enhancement factor z~ x = [1 — dReE^O)/^]. Especially, 
the former gives rise to the depairing effects of the Cooper pair due to the inelastic scattering, which would lower T c . 
For T <C T c , however, such a suppression does not occur because ImEo — > as T — ¥ 0. This property is a key for the 
colossal enhancement in H c2 for H || c-axis. 

We, next, consider the pairing interaction between the quasiparticles due to the strong spin fluctuations near the 
QCP. They are evaluated at the lowest order in g 2 Xa, 

V ss Mk,k') = -\g\(k-k') + \g 2 X {k + k'), (56) 
6 6 

V sS Mk,k') = ^9 2 x(k-k') + \g 2 x(k + k'), (57) 
6 3 

V sS , Ss {k,k') = -V srs Mk,-k r ), (58) 

and the other components are zero. These are directly derived from Eq.©. Although the applied fields might 
affect V, we neglect such an effect in this section. The 77-dependence of V can be included within our approach, 
if the field dependence of x(l) is clarified by some experiments. We also note that, in Eqs. ([55)) ^ ([551) . the spin-flip 
scattering processes are not included. They are expected to enhance the mixing of the spin singlet and the triplet 
superconductivity. These two neglected effects are discussed in Sec|Vl As noted in SecJITl the coupling constant g 
should be regarded as an effective one renormalized by the vertex corrections^^. 



A. H=0 case 



In this section, we study the gap function and the transition temperature at H = by solving Eq.(|40p. In this 
case, A does not depend on the center of mass R and Q is simplified as Q TT i{k) = G T (k)G T '(—k). Among the five 
irreducible representations of the point group Ci v for CeRhSi3 and CeIrSi3, the most stable symmetry of the gap 
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functions is Ai symmetry which is consistent with the previous study^ 3 - The ^-dependence of the singlet gap function 
is Agjngiet ~ cos(2/c z a), and that of the triplet gap function is A tr j p i ot ~ sin(fc 2 , J/ ), as will be discussed in detail in 
Sec IV AI In the previous study for H C 2 \\ z, we have neglected the triplet part of the gap function, because it is much 
smaller than the singlet one in amplitude. 2e In the present study, we take it into account and show that the results 
in the previous study are not changed. 

In Fig. ([2]), the transition temperatures for this Ai symmetric superconducting state for several 9 are shown as 
functions of g 2 Xo- T c saturates for large g 2 Xo because the strength of the pairing interaction and that of the depairing 




10 15 20 25 30 



FIG. 2: Transition temperatures T c /ti as functions of g 2 Xo/ti f° r several 8 at H = 0. 
0.002, 0.005, 0.01, 0.02 and 0.03 from the top to the bottom. 



The curves correspond to 9/ti 



effect through the normal selfenergy become comparable. Note that the dependence of T c on 9 is weak. 

The coupling constant g is fixed so that the calculated T c is of the same order as the experimentally observed T c . 
In CeRhSi3 and CeIrSi3, Kondo temperature is Tk ~ 50-100 (K)^ 2 - and the resistivity saturates around 200^300 
(K) ) 15 ' 17 which implies that the hopping integral t\ in our model is t\ ~ 50-100 (K). On the other hand, observed 
T c is T c ~ 1(K), that is, T c ~ 0.01ti-0.02i x . To reproduce this T c in the calculation, we fix g 2 Xo/h = 10-15. For 
these values, the system is in a strong coupling region. The renormalization factor averaged on the Fermi surface is 
z^ 1 1.7 for g 2 Xo/ti = 10, and z is not sensitive to 9, which is characteristic of the 3D AF spin fluctuations ! 21 ' 22 ' 57 ' 58 
Below, we mainly study the case of g 2 Xo/ti = 10 for which T c for the minimum 9 = 0.002ii is T c = 0.0139ii. Setting 
T c = 1.3(K), which is an averag ed value of T c for CeRh Si 3 ^ and CeIrSi 3 ,^ we have t x = 93.8(K). We also consider 
the case of g 2 Xo/ti = 15 for in-plane fields in Sec lIVCI In this case, similarly, we have T c — 0.0199ti = 1.3 (K) and 
ti=65.3(K). 



H II c-axis case 



In this section, we calculate the upper critical fields for H _L z. In this case, the parameter A which characterizes 
the anisotropy in the ii^i^-plane is A = 1. The other parameter Q which should be optimized is Q = because, for 
Q to be finite, the interband pairing on the split Fermi surface is required. However, such a pairing is energetically 
unfavorable. 

To study H C 2, we fix the strength of the coupling constant as g 2 Xo — lOij.. In this calculation, the admixture of the 
singlet and the triplet components of the gap functions is fully taken into account, which is neglected in the previous 
paper^i The results are almost unchanged from the previous ones even if we include the effect of the admixture. 
In Fig 121 H c2 || z curves as functions of T for several 9 are shown. The Pauli limiting field Hp is large because the 
Rashba SO interaction is strong, a = 0.5ii > T C (H = 0) ~ O.Olii. In such a case, the quasiparticles are easily paired 
on the same band under the applied field /ibH <C a. This holds generally and does not depend on the symmetry 
and the dominant parity of the gap functions for Rashba superconductors ^ 1 1 ' 1 2 The upper critical field is, therefore, 
mainly determined by the orbital limiting field i? rb- However, as seen in FigJ31 the orbital limiting field i/ rb is 
different from H C 2 calculated with both the Pauli and the orbital depairing effect being taken into account, especially 
for large H. It would be natural to think that this difference is a numerical artifact due to our choice of parameters. 
The magnitude of a used in the above calculations is not sufficiently large for high H regions. If one uses the large 



value of a/T c (0) , where T c w = ti exp[— l/(po5 Xo)] with the density of states at the Fermi level po, this difference 
may disappear. In fact, in the experimental data of H C 2 both in CeRhSi3 and CeIrSi3, no clear Pauli depairing effect 
can be seen, which implies that the Zeeman effect is effectively negligible in the compounds. However, to carry out 
the numerical calculations of H C 2 for the larger a/T c , we need the large size of the /c-mesh and a large number of 
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FIG. 3: H c2 (T) at g 2 Xo/ti = 10 for several 9. The dotted curve with triangles is Pauli limiting field H P for 9/ti = 0.002 
and the solid curves with circles are orbital limiting fields H OI h for 9/ti = 0.002,0.03. The dotted curves with squares are 
H C 2 curves including both the Pauli and the orbital depairing effects for 9/ti = 0.002, 0.005, 0.01, 0.02, 0.03 from the up to the 
bottom. The TZorb curve for 9/ti = 0.03 coincides with the H C 2 curve with squares. 

Matsubara frequencies. We have also calculated H c2 || z for g 2 Xo/ti — 15- Although, in this case, the Zeeman effect 
much affects H C 2 compared with the g 2 Xo/ti — 10 case, the qualitative behavior of _ff or b is unchanged. 

The calculated H C 2 show (i) strong 8 dependence and (ii) upward curvatures, and (hi) they reach ~ 30(T). All these 
characteristic behaviors well explain the experimental observations in CeRhSi3 and CeIrSi3 discussed in SecOJ 24 ' 25 
The physical reason for these characteristic behaviors in H C 2 \\ z is quite simple. Because, for H || z, H C 2 is 
determined mainly by the orbital depairing effect and the orbital limiting held H m h can be strongly enhanced by 
the spin fluctuations near the QCP. In the quantum critical regime, as T is decreased below T C (H — 0), the pairing 
interaction V oc £, 2 {T) is increased in magnitude while the inelastic scattering between electrons is suppressed and the 
quasiparticle damping is decreased, 7(T) = — ImS^(T) — ► 0. This contrasting behaviors of the pairing interaction and 
the depairing effect lead to the large enhancement in H C 2 for T — > near the QCP. On the other hand, as discussed in 
the next section, the Pauli limiting field Hp is not so strongly enhanced by the spin fluctuations at low temperatures. 
This is a key to resolve the apparent contradiction that although there are many heavy fermion compounds which are 
considered to be located near magnetic QCPs, they do not show such a huge H C 2 as in CeRhSi3 and CeIrSi3. In usual 
centrosymmctric heavy fermion superconductors, H C 2 is considered to be mainly determined by the Pauli depairing 
effect. Therefore, even if the system is close to the QCP, H C 2 is not anomalously enhanced. 

The pressure(#) dependence of H C 2(T — > 0) shows a remarkable feature as a result of above mentioned mechanism. 
We define normalized T c and H c2 as functions of 9, t c (9) = T C (H = 0,9)/T c (H = 0,9 = 9 U ) and h c2 {6) = H c2 (T = 
T m ,9)/H c2 {T = T m ,9 = 9 M ) where 9 M = 0.03ii and T m = 0.002ti. The normalized orbital limiting field h OT b is 
also defined in the same way. In Fig|2J t c ,h C 2 and /i or b are shown for g 2 Xo = 10ti . For h C 2, the dotted curve 
with triangles is calculated from i? rb, and the dotted curve with squares includes both the Pauli and the orbital 
depairing effect. The 6 dependence of t c is moderate, while those of both h C 2 and ft. or b are significant. As explained 
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FIG. 4: t c and h C 2 as functions of 9. The dotted curve with triangles is calculated from -ff or b, and the dotted curve with squares 
includes both the Pauli depairing and the orbital effects. The definitions of t c (9) and h C 2(9) are given in the text. 

above, these behaviors are understood as a result of the strongly enhanced pairing interaction and the suppression 
of the depairing effect at low temperatures in the vicinity of the QCP(6' = 0). Since, in CeRhSi3 and CeIrSi3, the 
SO interaction makes the superconductivity orbital limited, the huge H C 2 is a result of the interplay of the Rashba 
SO interaction and the electron correlations. Generally, such strong enhancement in the pairing interaction and the 
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suppression of the quasiparticle damping at low temperatures are crucial for orbital limited superconductors, because 
i/ or b is largely affected by the electron correlations compared with Hp. Therefore, the enhanced upper critical field 
can be considered as a universal property of the orbital limited superconductors near QCPs. This would be related 
to the recent experiments of H c2 || a-axis in UCoGe in which the relation between the superconductivity and the 
ferromagnetism has been discussed^ The observed H® 2 is huge ~ 15 (T) while T c ~ 1 (K)i 60 ' 61 This issue is now 
under investigation. 



C. H _L c-axis case 



We also study H C 2 for the case of H _L z within the same framework. Since in this case, the Fermi surface 
is distorted asymmetrically by the in-plane filed through the Rashba SO interaction, the Pauli depairing effect is 
significant, which implies that the higher Landau levels become important. Furthermore, the optimization parameter 
Q and A are nontrivial for H _L z. First, at a fixed H, we optimize A which corresponds to the anisotropy in the 
quasiparticle velocity of the two directions perpendicular to the applied field, or the anisotropy in the superconducting 
coherence length. Since A is characterized by the shape of the Fermi surface, the field dependence of A is very weak. 
We fix the optimized A, and then, optimize Q to have the maximum H C 2 for given temperatures. The optimal A is 
A -2.3. 

In FiglSJ we show H C 2 at g 2 xo = 10fi for two values of 9, 9yi/t\ = 0.03 and m /ti — 0.002. Each H c2 curve is 
calculated with a single Landau function for N = 0,1,2, respectively. True H c2 curve should be calculated by a 
superposition of the Landau functions. We have computed a H c2 curve by using a superposition of N = and N = 1 
Landau functions, and found that it almost coincides with the H c2 curve calculated by N = Landau function only. 
Therefore, H c2 is mainly determined by the N = Landau level, and the shapes of N = H C 2 curves for 9 = 9 m and 
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FIG. 5: H C 2 .1 z at g 2 xo/ti = 10 for 8 — 8 m (square symbols) and 8 = 8m (circle symbols). For each 8, three curves correspond 
to TV = 0, 1, 2 Landau levels from the top to the bottom. 

9m are similar. This pressure insensitivity is due to the weak dependence of the Pauli limiting field Hp on the electron 
correlations compared with H m ^. The ratio of the calculated value of H C 2(T — >• 0) _L z to that of H C 2(T — > 0) || z in 

the previous section is H^/H^ ~ 1/3 for 9 = 9 m . These behaviors in H c2 are consistent with the experiments . 24 ' 25 

We turn to the discussion of the modulation vector Q. Under the field [IbH <C a, the dispersion is changed as 
e T (k + q; H) ~ s T {k) +v T (k) ■ q + TfiB^(k) ■ H. In this situation, the momentum pair {kp T -\-Q T , —kp T ) is energetically 
degenerate on the one band, where kp T is the Fermi momentum for r-band and Q T satisfies v T ■ Q T = —t[1bC> • H. 
Note that the center of mass momenta of the pairs on each band satisfy Q + ~ — Q_, and the electrons on each 
band favor the center of mass momenta with opposite directions. Therefore, it is expected that for sufficiently strong 
H, each band favors each Q and the resulting superconducting state would be the Fulde-Ferrell-Ovchinnikov-Larkin 
(FFLO) stat o 35 ' 36 with A ~ r/ + exp[iQ + ■ R] + rj- exp[iQ_ ■ R]. For small H, however, the helical vortex state with 
A ~ r)exp[iQ ■ R] is considered to be stabilized in general noncentrosymmetric superconductors. So, the situation is 
different between the case of small H and that of large H. We discuss the f/-dependence of the modulation of the gap 
function qualitatively. When the SO split inter-band pairing which is small for a T c is neglected, the Eliashberg 
equation for the diagonal element of the gap A r is of the form, 



A T {a,H,Q) ~ p + (a,H)F r+ (a,H,Q,A+)+ P -(a,H)F T -(a,H,Q,A-) 

= \(P++ P-)(^r+ + Ft-) + \{P+ - P-)(J r T+ - Fr-), 



(59) 
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where p T is the density of states at the Fermi level for the r-band and T T ± is a function depending on (a, H, Q, A r ). 
The first term in Eg. (|5§1) is proportional to just the sum (p + + p_) and therefore, the electrons on each band 
contribute independently. In contrast, in the second term, the difference between the two bands plays important 
roles. The term is related to the magnetoelectric effect in the superconducting state due to the anisotropic SO 
interaction, which depends on the difference in the densities of states of the two bands. This effect is incorporated 
into the Ginzburg-Landau free energy as / mc oc H^K,^ v [ip* D v ^p + ip(D u ip)*] where = — i2eA fl and K,^ is the 
coefficient of the magnetoelectric effect ^^^S ~ 12 ' 27 i 34 In general noncentrosymmetric superconductors, / mo leads to a 
spatially modulated gap function with the modulation vector Q™ e oc H^K,^. In this state, the Cooper pair is formed 
by the states with kp T + Q m ° and —kp T momenta, not the (kp T + Q T , —kp T ) momenta. This effect arises even 
under very weak H . However, it is pointed out that in 3D Rashba superconductors in which only TC xy = —JC yx are 
nonzero, the phase exp[iQ mc ■ R] is absorbed into the Landau function as a spatial shift <fi(R) — > <fi(R — Rq) with a 
-H-dependent vector Rq£2.~M. Therefore, Q mo does not appear in physical observables like H C 2, although K, xy itself is 
nonzero. On the other hand, under a high field, the first term in Eq. (|59|) plays important roles for the optimization 
of Q. This situation is similar to that of the FFLO state in usual centrosymmetric superconductors. Since, as noted 
above, the first term of Eq. (|59|) is a sum of the independent contributions from two bands, it merely favors Q + or 
Q . Therefore, in the high field region, the candidate for the modulation vector is Q + and Q . If H is applied from 
zero to some large value for general noncentrosymmetric superconductors, we would see a continuous change of Q, 
from Q mc to Q ± . The threshold value H — H* at which Q changes from Q mc to Q ± depends on the details of the 
system. It is pointed out that H* becomes large as the orbital depairing effect increases. 

In our study, Q is determined so that H C 2 becomes maximum for a given parameter. We find that the optimized 
Q vector is parallel to a-axis, Q — (— Q,0,0) for H = (0, if, 0). In FigO the optimized Q along the H c2 curve for 
9 = 9 m is shown. We have two regions; the region where Q is quite small and the other with large Q. For the small Q 
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FIG. 6: The modulation Q along the H C 2 curve for g 2 Xo/ti = 10 and 6 = # m . In the low field region, the optimization of Q 
would be due to a numerical artifact and Q is shown by the dotted curve in this region. H* is the threshold value. 

region, although we have a systematic change of H C 2 with respect to Q and can optimize it, this dependence of H C 2 
on Q would be a numerical artifact because the change in H C 2 due to nonzero Q is infinitesimally small. Actually, the 
small Q region corresponding to the helical vortex state is spurious in a 3D Rashba superconductor because of the 
reason mentioned above. It is expected that the character of the stable vortex state is nothing but the character of the 
conventional vortex state with Q = in the region < H < H*. In contrast, for large H > H*, we have a finite Q. 
This large Q state would not be a direct result of the lack of the inversion center. Rather, it is stabilized by the pairing 
of the momentum (kp T + Q T , ~kp T ) electrons on each band. However, the contribution to the gap function from the 
second term in Eq. (|59")) is not negligible, resulting in the shift of the degeneracy between Q + and Q-. Therefore, we 
expect that, in this high field region, the FFLO state with the gap function A ~ rj + cxp[iQ + • R] + r\- cxp[iQ_ • R] 
can be realized. We have performed the calculations for other parameters, and confirmed that the threshold H = H* 
for the two region depends on the effective mass m e g. As m e g becomes larger, H* decreases, and vice versa, which 
means that the orbital depairing effect plays important roles for the determination of H*. To discuss the stability 
of such a state, we need to compute the free energy in the superconducting state, and it is beyond the linearized 
calculation of H C 2 performed in the present study. 

As mentioned in Sec lIIIl g 2 Xo cannot be determined uniquely in our theory, and we can change the value of g 2 Xo 
within the range for which the value of T c is consistent with the experiments. In the following, we study the case of 
g 2 Xo — 15ii which gives T C (H — 0, (9 m ) = 0.0199£i. In this case, the value of t\ is t\ — 65.3(K) and the effective mass 
of the cyclotron motion m e s is large compared with the t\ = 93.8(K) case. We show H C 2 for 9 = 9 m , 9m in Fig{7j The 
left panel shows H C 2 curves calculated with the single Landau functions for N = 0,1,2, and the right panel shows 
H C 2 curves calculated with the superpositions of the N = and the N = 1 Landau functions. For 9 — 9 m , the H C 2 
curve calcuated with the N = 1 Landau level is larger than that with the N = Landau level at a low temperature 
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FIG. 7: H C 2 1 2 at g 2 Xo/ti = 15 for = # m (square symbols) and = 0m (circle symbols). Left panel; For each 0, three 
curves correspond to N = 0, 1,2 Landau levels from the top to the bottom. Right panel; H C 2 calculated with the use of the 
superposition of the N = and N — 1 Landau functions. 



region. In such a region, higher Landau levels become important, and the gap function can have the nodal structure 
in real space due to the nodes of the higher Landau functions ! 32 ' 33 The H C 2 curves calculated with the use of the 
superpositions of the N = and the N = 1 Landau levels almost coincide with the N = H c i curve for low H 
and the N = 1 H C 2 curve for high H, respectively. In the case for 6 — m , the higher Landau levels becomes more 
important than the case for = 8m, because the orbital depairing effect is largely suppressed and the electrons are 
strongly paired near the QCP. 



V. SPIN-FLIP SCATTERINGS AND FIELD DEPENDENCE OF SPIN FLUCTUATIONS 



In the calculation shown in the Sec lIVl we have neglected two important effects, the spin-flip scattering processes 
in the pairing interaction and the field dependence of the spin fluctuations. Regarding the former, in the noncen- 
trosymmetric systems, there always exist spin-flip scattering processes which are not included in Eqs. (|5o]) ^ (f55)) . It 
was pointed out that they can enhance the mixing of the singlet and the triplet superconductivity; 37 ' 38 and also, the 
effective strength of the Pauli depairing effect depends on the ratio of the admixture of the gaps for H _L zM- Another 
important point which is neglected in the calculation in Sec lIVl is the field dependence of the susceptibility. Because 
the observed H C 2 is over 20(T) for c-axis in CeRhSi3 and CeIrSi3, one might think that the spin fluctuations are 
suppressed by such a strong magnetic field, although we have assumed in Eqs. (l56|) ^ (|58|) that the spin fluctuations are 
not strongly affected by the magnetic field. These two points are examined in this section, and it is concluded that 
the neglect of them is a legitimate approximation and the calculated results in Sec lIVl are qualitatively unchanged 
even if we take into account the two points. 



A. spin-flip scatterings in pairing interaction 



In this section, the effects of the spin-flip scattering processes in the pairing interaction on the superconductivity 
are examied. Through a spin-flip process, such as the scattering process in which spin f! particles are scattered as 
spin ft particles, the singlet and the triplet pairing states are mixed directly. It is pointed out by several authors that 
this effect can enhance the admixture of the parity even and odd pairing. 37,38 It is also discussed that, for in-plane 
fields, the effective strength of the Pauli depairing effect depends on the ratio of the triplet gap function to the singlet 
gap function^ 3 - In the following, we show that, in CeRhSi3 and CeIrSi3, the admixture of the gap functions is not so 
strong even if we include the spin-flip scattering processes in the pairing interaction. 

To investigate the effect of the spin-flip, we use the single band Hubbard model 

H = *y] c\e {k)c k + a^2 c{Co(k, H) ■ ac k + U ^ n iT n 4 . (60) 

k k i 

Here, as in eq.([T]), Cfc s is the annihilation operator of the Kramers doublet of the heavy electrons which are formed 
through the hybridization with the conduction electrons. The dispersion relation eo{k) and the Rashba SO interaction 
C (k,H) are defined in Eqs.Q and ©. We fix the parameters as {t x ,h,hM>n,a) = (1.0,0.5,0.3,0.025,0.975,0.2) 
in this section and the next section. The pairing interaction is evaluated by the random phase approximation (RPA) 

V a p !a '0'{k, k') = U a p tCt >p> + [U + \(k + k')U + ]/3> a> a > - [U-x(k - k')U-] aa ' t pp (61) 
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where the matrices are defined with the notation 

Mffff -Mfffi 

]\/T. ... i\/r. , . , l\/f 

M = 



Mj-|-|~|- Mf-j-fj, Mfftf -^ttU 

Mf4.f1 Mf4f4 Mf4_4f Mf444 

Mj,fff M^ffi Mj,f4,f M|f44 

4-4'TT "J'4''t'4' "4*4''4"'T — i — X — 



M 



(62) 



The matrices U, U+ and [/_ arc defined as 
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The susceptibility x(q) within RPA is 

X(q) = X°(q)[l-U + X°(q)} 



(63) 



(64) 



(65) 



(66) 
(67) 



Equation tfBTj) includes the spin-flip scattering processes, even for V SSS s and V S s S s as the virtual scattering processes. 
The matrix interation V is characterized by the susceptibility and, in the limit of a — > 0, it coincides with 
Eqs. (f5"6l) ^ (l5"8|) if we neglect the onsite repulsive term U and the charge susceptibility terms. As shown in Sec lVBl 
x(q) has a peak around q ~ (±0.57r, 0, 0.57r) and q ~ (0, ±0.57r, 0.57r) which is consistent with the neutron scattering 
experiments for CeRhSi3, and the q-dependence of x(q) is almost the same as the phenomenological x(q) defined by 
Eq.©. 

To discuss the effect of the spin-flip processes on the admixture of the singlet and the triplet gap functions, we solve 
the Eliashberg equation within the weak coupling approximation, 



A aa , (ty = -—^2 Vaa'PP' ( k i k')gp 7 p> 7 > (k') A lY (k r ) 



(68) 



where Vaa'pp' k 1 ) is calculated with the use of Xaa< 0p>{iVn = 0, q) and g^ 7 ^' 7 '(fc) = T^^^ G° 7 (fc)G°, 7 ,(— fc). 
The pairing interaction V consists of two parts, V CO n corresponding to the spin conserving scattering processes, and 
Vflip including the spin-flip scattering processes. For convenience, we introduce the 4-component <i-vector of the gap 
function, using the identity matrix ctq and the Pauli matrices <r, 



(69) 



where do and d are, respectively, the singlet part and the triplet part of the gap functions. We calculate d M (fc) for two 
cases; (i) one where all the elements of the interaction matrix V aa i /3/3>(k, k') are fully taken into account, and (ii) the 
other where the terms V con including the spin-flip processes are neglected. We fix the parameters as U — 3.5ii,T = 
0.02ti. For these parameters, the eigenvalues of the Eliashberg equation are 0.95 ~ 1.05. The gap functions for case 
(i) are shown in Fig(5]for the singlet gap function do and in Fig |H] for the triplet gap function d\. For case (ii), FigfTUl 
and FigfTTI show do and d\ , respectively. In both cases, the singlet gap function is approximately do(k) ~ cos(2k z a). 
The ratio of the triplet gap function to the singlet gap function defined as r = |max{di(fc)}|/|max{<io(fc)}| is about 
r ~ 1/10 for case (i) and r ~ 1/30 for case (ii). Although r is enhanced by the spin-flip processes, it still remains 
small in our system. We have performed similar calculations for various Fermi surfaces, and found that, generally, 
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E3 



k x =0 




FIG. 8: The singlet gap functions dd(fe) in k z — it/2 plane (left) and k x = plane (right). The spin-flip scattering processes 
are fully taken into account. The broken curves represent the Fermi surface. 



k z =7t/2 



k x =0 




FIG. 9: The triplet gap functions di(fe) in k z — tt/2 plane (left) and k x — plane (right). The spin-flip scattering processes 
are fully taken into account. The broken curves represent the Fermi surface. 



the spin-flip processes can enhance r. However, the value of r depends on the details of the system. In CeRhSi3 and 
CeIrSi3, we conclude that the admixture of the gap functions is small. Therefore, the results in Sec II VI where we have 
neglected the spin-flip processes in the pairing interaction are supported. Finally, we note that the fc-dependence of 
the triplet gap function is largely affected by the spin-flip processes. If we write V = V CO n + cVfli p with a tuning 
parameter c, the change in with respect to c is continuous, although d\ for case (i) and (ii) look quite different 
from each other. This difference in the fc-dependence in d\ is not important for the discussion of H C 2 in Sec lIVI 



k z =n/2 



k x =0 




FIG. 10: The singlet gap functions do(k) in k z = tt/2 plane (left) and k x — plane (right). The spin-flip scattering processes 
are not taken into account. The broken curves represent the Fermi surface. 
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k z =7t/2 



k x =0 




FIG. 11: The triplet gap functions di(fc) in k z = n/2 plane (left) and k x = plane (right). The spin-flip scattering processes 
are not taken into account. The broken curves represent the Fermi surface. 



B. field dependence of spin fluctuations 

The field dependence of the spin fluctuations is another effect which has been neglected in Sec lIVI Experimentally 
observed H C 2 is so large especially for H C 2 \\ z that one might think that the susceptibility x(q) is affected by the 
applied field and the spin fluctuations are weakened. We show, however, that the effect of the applied field on x(q) is 
strongly suppressed by the Rashba SO interaction. This is because the Rashba SO coupling tends to fix the direction 
of the spins on the Fermi surface depending on k- vectors, which competes with the Zeeman effect. As a result, the spin 
fluctuations in CeRhSi3 and CeIrSi3 are robust against the applied magnetic field up to the strength of the Rashba 
SO interaction, [IbU ex. 

We compute x{l) under finite fields H = (0, H y ,0) or (0, 0, H z ), 

XM = [^dre^ (T^(r)S^(O)) 

•'o (70) 
1 

where XaPP'a'il) 1S evaluated within RPA used in Sec IV Al We fix U and T, as in the previous section, U = 3.5ii 
and T — 0.02ii. In Fig|12j iJ-dependence of Xxx,Xyy an d Xzz is shown for H = (0,0, H z ). At H = 0, Xxx{0,Q x ) = 



q 2 =ii/2 




(0,0) (it,0) (0,71) (jc,j[) (0,0) 



q 2 =ii/2 




(0,0) (i,0) (0,n) (j[,x) (0,0) 



q z =Jt/2 




(0,0) (n,0) (0,Jt) (j[,jt) (0,0) 



FIG. 12: Perpendicular filed h z = fiBH z /ti dependence of the susceptibility XtwiWn ~ 0, q) at q z = tt/2 for U = 3.5t\ and 
T — 0.02ti. The strength of the Rashba SO interaction is a = 0.2ti. Xxx(\eft), Xyy( cen ^ eT )^ an d Xzz( r ight,)- 

Xyy(®- Qy) > Xzz{0, Q x y ) are satisfied, which is consistent with the result of the neutron scattering experiments in 
CeRhSi3 that the antiferromagnetic moment is in afr-plane.— Here, Q x ~ (±0.57r, 0, 0.57r) and Q y ~ (0, ±0.57r, 0.57t). 
Note that Xxx(0, Q x ) > Xyy(0, Q x ) is satisfied because of the spin-flip scattering processes. However, the anisotropy in 
the non- interacting x° M (g) is of the order of a/sp <C 1, and therefore, the anisotropy in Xnnio) including the electron 
correlation effect remains irrelevant for the discussion of H C 2 even near the QCP. As mentioned above, for fisH < a, 
{Xfift} are almost unchanged. For [IbH > a, only Xzz is suppressed. The robustness of for [IbU ^ a is a 

general feature of the noncentrosymmetric systems, since the spins for every fc-point are fixed by the anisotropic SO 
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interaction in that region. These calculated results support the legitimacy of our neglecting the field dependence of 
the pairing interaction for the calculation of H c2 . Although there is no direct observation of the strength of the SO 
interaction, it is expected to be pretty large, a > hbH c2 ~ 30//b(K). Therefore, in CeRhSi3 and CeIrSi3, the spin 
fluctuations remain so strong under applied magnetic fields that H C 2 is strongly enhanced. 

The same robustness also exists for H = (0, H y , 0) for which the Fermi surface is distorted anisotropically. FigurdTSl 
shows the field dependence of Xw All of Xnn are almost unchanged for [isH < a, and \yy is suppressed for [ibH > a. 
These behaviors are basically the same as those for H || z, and assert the robustness of the spin fluctuations against 




(0,0) (it,0) (0,ii) (it,Jt) (0,0) (0,0) (i,0) (0,n) (it,Jt) (0,0) (0,0) (11,0) (0,31) (Jt,jt) (0,0) 

q q q 



FIG. 13: In-plane filed h y = /isHy/ti dependence of the susceptibility x^{0,q) at k z = tv/2 for U = 3.5ti and T = 0.02ti. 
The strength of the Rashba SO interaction is a = 0.2ti. Xxxfleft), \w (center), and Xzz (right). 

the in-plane field. 

From the above results, we see that neglecting the field dependence in the pairing interaction for any field direction 
is legitimate provided [LbH C 2 ^5 a - Although H C 2 is huge in CeRhSi3 and CeIrSi3 especially for c-axis, the condition 
HbH C 2 < a is expected to be satisfied. Therefore, the discussion of H C 2 in Sec|lV]is not changed even if we consider 
the field dependence of the spin fluctuations. 



VI. SUMMARY 



We have discussed the normal and the superconducting properties in noncentrosymmetric heavy fermion super- 
conductors CeRhSi3 and CeIrSi3. We have shown that the T- linear dependence of the resistivity above T c observed 
experimentally is naturally understood within the 3D spin fluctuations near the AF QCP. 

For the superconducting state, we have derived a formula from the Eliashberg equation in real space. The formula 
enables us to treat the Pauli and the orbital depairing effects on an equal footing. Furthermore, by using it, we can 
calculate H C 2 for strong coupling superconductors with general Fermi surfaces. We have calculated H C 2 with the 
formula and have well explained the observed features of H c2 in CeRhSi3 and CelrSi3. For H \\ z, Hp is infinitely 
large due to the Rashba SO interaction, and H c2 is determined by -fforb- As temperature is lowered and the system 
approaches the QCP, the pairing interaction becomes larger while the quasiparticle life time becomes longer, which 
results in the huge -fforb — H c2 with the strong pressure dependence. The enhancement of the orbital limiting field 
near QCPs by this mechanism would be universal. We have also discussed the case for H _L z. In this case, the 
Pauli depairing effect is significant because of the asymmetric distortion of the Fermi surface, and resulting H C 2 is 
moderate against the pressure. The FFLO state can be stabilized for a large H region, although such a region is very 
small. The features of the calculated H C 2 for both H || z and H _L z are in good agreement with the experiments. 
This consistency supports the scenario that the superconductivity in CeRhSi3 and CeIrSi3 is mediated by the spin 
fluctuations near the AF QCP. 

In the last section, we have checked the legitimacy of our approximation used for the calculation of H C 2 . In CeRhSi3 
and CeIrSi3, the admixture of the singlet and the triplet gap functions are small even if we take into account the 
spin-flip scattering processes in the pairing interaction. In noncentrosymmetric systems, the spin susceptibility is 
robust against the applied magnetic fields /is-ff < a. For this reason, the spin fluctuations near the AF QCP in 
CeRhSi3 and CeIrSi3 remain strong even under a large magnetic field ~ 30 (T). Therefore, the above mentioned 
results for H C 2 is not changed if we refine our approximation used in the calculation of H C 2- 
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